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Abstract 

Background: Little is known about the spatiotemporal pattern of bacillary dysentery (BD) in China. This study assessed the 
geographic distribution and seasonality of BD in China over the past two decades. 

Methods: Data on monthly BD cases in 31 provinces of China from January 1990 to December 2009 obtained from Chinese 
Center for Disease Control and Prevention, and data on demographic and geographic factors, as well as climatic factors, 
were compiled. The spatial distributions of BD in the four periods across different provinces were mapped, and heat maps 
were created to present the seasonality of BD by geography. A cosinor function combined with Poisson regression was used 
to quantify the seasonal parameters of BD, and a regression analysis was conducted to identify the potential drivers of 
morbidity and seasonality of BD. 

Results: Although most regions of China have experienced considerable declines in BD morbidity over the past two 
decades, Beijing and Ningxia still had high BD morbidity in 2009. BD morbidity decreased more slowly in North-west China 
than other regions. BD in China mainly peaked from July to September, with heterogeneity in peak time between regions. 
Relative humidity was associated with BD morbidity and peak time, and latitude was the major predictor of BD amplitude. 

Conclusions: The transmission of BD was heterogeneous in China. Improved sanitation and hygiene in North-west China, 
and better access to clean water and food in the big floating population in some metropolises could be the focus of future 
preventive interventions against BD. BD control efforts should put more emphasis on those dry areas in summer. 
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Introduction 

Diarrhoea is an important cause of disease burden worldwide, 
particularly for resource-limited countries and for children under 
age of five [1,2,3,4]. Bacillary dysentery (BD), a severe form of 
shigellosis, accounts for approximately five percent of diarrhoeal 
episodes [5,6], and it contributes to over a million deaths annually 
worldwide [6]. Due to the increasing resistance of shigella to 
multiple antimicrobials [7], it is essential to develop effective 
strategies to control and prevent BD. China has made rapid 
economic progress in the past decades, and the morbidity and 
mortality of BD have been declining progressively [8], partially 
due to the substantially improved water supplies, sanitation and 
hygiene [9]. Nevertheless, it still shoulders a considerable burden 
of disease from BD [5]. 

Recognition of the geographic distribution of BD can guide 
future health resource allocation, and exploring BD seasonality 
and understanding the season-specific risk factors can shed new 
light on the development of future vaccination programs. So far, 
few data on the spatial pattern of BD in China are available, and 



existing literature offers limited information on the pattern and 
potential drivers of BD seasonality in the whole China. This study 
used data on the monthly BD cases in 3 1 provinces of China from 
1990-2009 and attempted to: I) explore the geographic distribu- 
tion of BD in China; II) characterize the seasonality of BD by 
geography in China; and III) identify the putative drivers of BD 
morbidity and seasonality. 

Materials and Methods 

Data collection 

BD is a legally notifiable infectious disease in China. All clinical 
and hospital doctors are required to report cases of bacillary 
dysentery to local Center for Disease Control and Prevention 
(CDC). From 1990 to 2003, local CDCs reported the cases to 
China CDC once every month. From 2004 to 2009, the reporting 
system became internet based and doctors reported BD cases 
through Chinese Information System for Diseases Control and 
Prevention (CISDCP) [10]. Monthly data on BD cases in China 
from I s January 1990 through 31 s December 2009, nationally 
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Figure 1. The decomposed distribution of monthly bacillary dysentery morbidity in China, from January 1990 to December 2009. 

doi:1 0.1 371 /journal.pone.01 04329.g001 



and provincially, were obtained from China CDC. The onset 
month and place of the cases were included in the data. Data on 
BD included both probable and lab-confirmed cases as defined by 
the national BD standard. Illness is characterized by tenesmus 
symptoms, pus stool, mucus stool, watery stool, or loose stool. 
Laboratory criteria for diagnosis are based on the positive isolation 
of Shigella bacteria. 

Data on the demographic (population and per capita gross 
regional product (PGRP)) and geographic (latitude and longitude) 
characteristics and climate variables (monthly average mean 
temperature, relative humidity and rainfall) in different provinces 



were retrieved from China Statistical Yearbooks [11]. Ethical 
approval was obtained from the Human Research Ethics 
Committee of Queensland University of Technology (Australia) 
and the Research Institutional Review Board of Public Health of 
Shandong University (China) prior to the data being collected. 
Patient information was de-identified and thus no written 
informed consent was obtained. 

Statistical analysis 

A seasonal decomposition analysis was conducted to assess 
whether there was a long-term trend and a distinct seasonality of 



Table 1. Summary statistics for monthly incidence of bacillary dysentery (1/10 000) in China, from 1990-2009. 





Time 


Mean 


Standard deviation 


Minimum 


Maximum 


Morbidity (1/10 000) 


Period 1 (1990-1994) 


0.736 


0.654 


0.100 


2.667 


Period 2 (1995-1999) 


0.493 


0.370 


0.099 


1.417 


Period 3 (2000-2004) 


0.318 


0.207 


0.095 


0.704 


Period 4 (2005-2009) 


0.230 


0.152 


0.070 


0.652 



doi:1 0.1 371 /journal.pone.01 04329.t001 
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Figure 2. The spatial patterns of bacillary dysentery morbidity across four periods in China, from 1990-2009. 

doi:1 0.1 371 /journal.pone.01 04329.g002 



BD in China [12]. Descriptive statistics of BD morbidity in China 
in four periods (1990-1994, 1995-1999, 2000-2004 and 2005- 
2009) were analysed to describe the dynamics of the disease. The 
spatial distributions of BD in the four periods across different 
provinces were also mapped. 

In terms of BD seasonality, we created heat maps to present the 
peak and trough times of BD in each province [13,14]. A cosinor 
function with Poisson regression was used to quantify the peak 
time, trough time and annual amplitude of BD [15]. For a better 
comparison between different provinces, the annual amplitude was 
expressed as a proportion of mean case counts. Regression tree 
analysis was conducted to identify the potential drivers of 
morbidity and seasonality of BD [16]. Visual maps were created 
using ArcGIS 9.3 (ESRI Inc., Redlands, CA, USA), and aU other 
analyses were performed using R (version 2.15), with the "season" 



package (version 0.3-3) to conduct the cosinor analysis and 
"rpart" package to conduct the regression tree analysis. 

Results 

The summary statistics for BD morbidity in China over the four 
periods, including mean, maximum and minimum values, were 
presented in Table 1 . The average values of BD morbidity over 
the four periods were 0.736, 0.493, 0.318 and 0.230 (per 10 000), 
respectively. The maximum and minimum values of BD morbidity 
also experienced a gradual decrease across the four periods. 
Figure 1 presented the 20-year decomposed distribution of BD 
cases, revealing a decreasing long-term trend and a distinct 
seasonality. 

Figure 2 showed the spatial patterns of BD morbidity over the 
four periods, illustrating that BD morbidity experienced a sharp 
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Figure 3. Heat maps of bacillary dysentery epidemiology data in the whole China (A) and 31 provinces of China (B), from 
1990 to 2009. 
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decrease in most regions during the 20 years. In North China 
(Beijing, Tianjin, Hebei, Shared and Inner Mongolia), a significant 
decrease in BD morbidity was observed. In North-east China 
(Liaoning, Jilin, and Heilongjiang), East China Jiangsu, Zhejiang, 
Anhui, Fujian, Jiangxi, and Shandong) and Central and South 
China (Henan, Hubei, Hunan, Guangdong, Guangxi and 
Hainan), gradual declines in BD morbidity occurred in the 20 
years. In South-west China (Sichuan, Guizhou, Yunnan, Chongq- 
ing, Tibet), we found BD morbidity declined sharply (from very 
high to very low) in Sichuan and Yunnan from the first period 
(1990-1994) to the second period (1995-1999), while BD 
morbidity in Tibet and Guizhou remained quite high in the first 
two periods (1990-1999). BD morbidity in most provinces of 
North-west China (Shaanxi, Gansu, Qinghai and Xinjiang) did 
not decrease much until the last period (2005-2009), and BD 
morbidity in Ningxia was consistently high during the 20 years. 

Figure 3A revealed a relatively constant seasonality in China 
from 1990-2009, indicating that BD was more likely to peak in 
summer (July-September), even though there were two exceptions 
(1990 and 2008). Figure 3B showed the seasonality of BD by 
geography (listed by latitude), revealing that BD appeared to peak 
from July-September at high latitudes, and peak from July- 
October at intermediate latitudes. At low latitudes, BD peak time 
varied gready between provinces (especially Yunnan and 
Guangxi), from May-October. Summer and autumn were the 
most pandemic seasons (Figure 4), and the morbidity difference 
between summer and autumn during 2000-2004 was higher than 
other three periods. Figure 5 presented the spatial pattern of BD 
amplitude, showing that higher latitudes (Northern China) were 
more likely to have greater amplitudes. 

The summary statistics for the demographic and geographic 
characteristics, climatic factors, and BD seasonal parameters of the 
3 1 provinces in China, which were included in the regression tree 
model, were presented in Table 2. The potential drivers of BD 



morbidity and seasonality were unveiled in Figure 6. Relative 
humidity and longitude were significandy associated with BD 
morbidity. Provinces with dry (monthly average relative humid- 
ity<58.5%) climate shouldered greater BD burdens. Relative 
humidity was also associated with BD peak time, with provinces 
with monthly average relative humidity <69.5% peaking earlier. 
PGRP was associated with BD trough time, and provinces with 
higher PGRP reached trough time earlier than those with lower 
PGRP. Latitude was the primary predictor of BD amplitude, and 
the two latitudinal thresholds were 29.85°N and 37.1°N 

Discussion 

This study examined the spatiotemporal pattern of BD in China 
from 1990 to 2009. China has experienced a significant decrease 
in BD morbidity over the past two decades, although the decrease 
was not uniform between regions and Beijing and Ningxia still 
shouldered heavy BD burdens in 2009. In North-west China, BD 
morbidity declined slowly, compared with East and Central 
China. The majority of BD cases occurred between July and 
October, even though heterogeneity in BD peak and trough times 
existed between regions. Relative humidity was found predictive of 
BD morbidity and peak time. PGRP helped distinguish BD trough 
time, and latitude was the major driver of BD amplitude in China. 

BD risk map in China has been shrinking in the past two 
decades, due largely to the fast progressing economy, much 
improved hygiene and better access to sanitary water and food [9] . 
However, BD still remains the first or second dominating 
infectious disease in many regions [17,18]. In this study, we found 
that BD morbidity was unacceptably high in Beijing and Ningxia 
throughout the whole study period (1990-2009), indicating that 
not only regions with low economic development had high BD 
infection. Previous studies indicated that the high BD morbidity in 
Beijing can be attributable to the relatively low under-reporting 
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Figure 4. The distributions of bacillary dysentery morbidity across four seasons in China, from 1990 to 2009 (1 = spring, 
2 = summer, 3 = autumn, 4 = winter). 

doi:10.1371/journal.pone.0104329.g004 
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Figure 5. The relative amplitude of bacillary dysentery 
seasonality in China. 
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rate [17], the poor street food sanitation in some districts [19], and 
the emergence of shigellosis strains resistant to multiple antimi- 
crobials resulted from overuse of antibiotics [17,19,20]. In China, 
over the past decades, millions of people in rural areas, especially 
physical workers (e.g., builders), migrate to metropolis like Beijing, 
working for a better salary. Majority of the floating population live 
in areas with poor access to clean water, which may also cause the 
high BD morbidity in Beijing [20]. Prior research offers limited 
information on BD in Ningxia. Wang et al. analysed BD in 
Ningxia from 1958 to 1997 and found that students and farmers, 
and people living in places with close proximity to river had 
increased risk of BD [18]. People living nearby rivers are more 
likely to use river water for bathing, washing, and cleaning of 
eating, which might expose them more to sewage compared with 
people who live far away from river [2 1] . 

In this study, we found BD predominantly peaked between July 
and October in China, corresponding to the BD peak time in 
Vietnam [22]. Geographical diversity in BD peak time was also 
found, with BD at high latitudes (northern China) mainly peaking 
between July and August, and BD at low latitudes (southern 
China) peaking between May to October. Existing science offers 
limited information on the drivers of special BD seasonality in 
Yunnan and Guangxi, and we speculated that it may be because of 
the different monsoon seasons in these regions [23,24]. The motive 
of this study was to assist BD surveillance system and help develop 
future early warning system. In light of the diversity of BD peak 
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Figure 6. Predictors of bacillary dysentery morbidity and seasonality in China. 
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time by geography, people living in northern China should be 
made aware of the high BD risk from July to August and BD 
surveillance in southern China should be extensively strengthened. 
Greater BD amplitude of seasonality at higher latitudes was also 
detected in this study, indicating that some latitude-related factors, 
such as the reduced sunlight and its possible effect on Vitamin D 
level [25], may contribute to the BD seasonality in China [26]. 

The regression tree approach we used in this study avoided the 
multi-collinearity among different meteorological factors (temper- 
ature, rainfall and relative humidity), and the results suggested that 
relative humidity may impact BD morbidity, echoing to prior 



studies examining the relationship between climate factors and BD 
[27,28]. The underlying reason may be because humidity affects 
BD pathogens [27] and normally high relative humidity comes 
with high rainfall which reduces contaminated food and water 
[28]. In terms of the putative predictors of BD seasonality, we 
found relative humidity and PGRP were associated with BD peak 
time and trough time, respectively, and latitude was the primary 
predictor of BD amplitude, unveiling that BD seasonality may 
result from the interaction between climate, socioeconomic factors 
(which is associated with population's access to clean water and 
food) and population immunity to specific pathogens. 
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This study is, to our best knowledge, the first study to look at 
spatial pattern and seasonality in the whole China. The risk areas 
and peak and trough times we identified will assist BD 
surveillance. Several limitations of this study should be acknowl- 
edged. First, under-reporting of BD morbidity might result in 
information bias to some extent due to the use of surveillance data. 
Second, the provincial-level data we used are not ideal in 
identifying the small-scale BD high risk areas. Third, access to 
health care may vary across provinces, and the difference in 
reporting that may exist across provinces, which may result in 
information bias. 

Conclusions 

China has undergone a substantial decrease in BD morbidity 
from 1990-2009, but the progress is heterogeneous between 



regions. BD in China peaked mainly between July and October, 
but the peak month varied by place and time. People in western 
China, especially those living in regions with low socioeconomic 
status, were more likely to be attacked by BD. BD control and 
prevention interventions should focus more on the regions with 
large floating populations or less-developed economy, and BD 
surveillance in southern China needs to be strengthened. 
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